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We study the impact of very Ught sterile neutrinos (Am^(,^,, £ [1,10] x 10~^ eV^, sin^ 20nGw < 
10~^) on upcoming ^13-driven reactor antineutrino experiments like Double-CHOOZ and Daya Bay. 
Oscillations driven by these vales of Arrinew affect data in the near and far detectors differently 
and hence potentially modify the capability of these experimental setups to constrain and measure 
sin^2Si3. We find that the hypothesis Sncw 7^ negatively impacts one's ability to either place 
an upper bound on sin'^ 26ri in the advent of no oscillation signal or measure sin^ 2613 if a ^13- 
driven signal is observed. The impact of sterile neutrino effects, however, depends significantly 
on one's ability to measure the recoil positron energy spectrum. If sin'^ 2Sn(,w ^ 10"^, upcoming 
QQ ' Si3-driven reactor antineutrino experiments should be able to measure sin^ 20new and Arrinew, along 

\ with sin^ 2^13, as long as one is sensitive to distortions in the recoil positron energy spectrum in the 

. near (and far) detectors. 

(N ; 

(1^ . I. INTRODUCTION 

c/^ : 

^\ ' Next-generation reactor antineutrino experiments are currently under construction (for details on concrete projects 
(N ; see Guided by the results of current neutrino oscillation experiments ^ , these are aimed at observing electron 

antineutrino disappearance driven by the "atmospheric" mass-squared difference, Am^g, and hence measuring the 
(— I \ still elusive 6*13 mixing angle (c/. Sec. HIl). 

O ^1 In order to significantly improve on previous reactor antineutrino experiments, next-generation experiments will 
make use of a two-detector setup. A far detector is to be placed an optimal distance away from the reactor core so 
as to maximize 0i3-driven electron antineutrino disappearance, while a near detector is placed close to the reactor 
^ ' core in order to measure the "unoscillated" electron antineutrino flux. Assuming differences between near and far 
detectors are understood at the few permille level and that enough statistics are accumulated, one ultimately aims at 
being sensitive to sin^ 26113 values around 1%. 
" >^ ' The setup summarized above relies on the assumption that all neutrino oscillation phenomena are properly described 
in terms of three massive neutrinos which interact under the well-known weak interactions. Current data reveal this to 
be an excellent assumption, but the possibility of other "subleading" effects remains. Here, we consider the possibility 
I that light sterile neutrinos mix slightly with active neutrinos and hence modify the standard neutrino oscillation 
10 ■ picture. 

, Sterile neutrinos — gauge singlet fermions — are among the simplest extensions of the standard model of particle 
physics. They may be an integral part of the physics responsible for neutrino masses 0, Q , and may be a component 
00 \ of the dark matter in the Universe Theoretically speaking, nothing is known about sterile neutrino masses, 
■ and very light sterile neutrinos (say m,y <^ 1 keV) are as natural (as defined by 'tHooft) as very heavy ones (say 
TOjy ^ 10^" GeV). Experiments, therefore, provide almost all the unbiased information we have regarding sterile 
neutrino masses. We briefly discuss experimental constraints in Sec. HIl 

We concentrate on sterile neutrinos that could qualitatively affect the interpretation of 0i3-driven reactor neutrino 
experiments. We find that this can happen for very light sterile neutrinos that introduce to the three-neutrino- 
oscillation picture a new mass-squared difference of order (1 — 10) x 10~^ eV^. For example, the "hypothetical 
presence" of such sterile states hinders the ability of these experiments to rule out certain values of 6*13 if no evidence 
for oscillations is observed. On the other hand, small sterile neutrino effects, if present, may lead to very significant 
oscillatory effects in the near detector even in the limit of vanishing ^13. Detailed results are presented in Sec. IIIIl 
along with a study of how well ^13 and sterile mixing parameters can be simultaneously measured with next-generation 
reactor neutrino data. 

Before proceeding, we'd like to highlight our goals and the limitations of our analysis. We are interested in discussing 
the fact that sterile neutrinos can significantly alter the interpretation of the comparison between near and far detector 
data in reactor antineutrino experiments. We would also like to point out that the near detector — because it is 
expected to collect hundreds of thousands of neutrino scattering events — may play an active role in revealing new, 
unexpected physics. With this in mind, our simulations, discussed in more detail in Sec. IIIIl are not aimed at 
realistically describing experimental setups or quantitatively gauging their capabilities. Qualitatively, however, we 
believe that all of the effects discussed in this manuscript will manifest themselves once real experimental data is 
analyzed. Our main message, along with our results, is summarized in Sec. IIVI 
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II. VERY LIGHT STERILE NEUTRINO EFFECTS AT REACTOR EXPERIMENTS 



We are adding to the three active neutrinos a fourth sterile state and are hence faced with four neutrino mass 
eigenstates. The electron neutrino can be expressed as a linear combination of f^, i = 1,2,3,4: v^, = ^Yl,i^ei^i- 
We order the neutrino masses as follows (for a detailed discussion see 0). 1^4 is the "mostly sterile" state (|t^s4p ^ 



|C/si,s2,s3p) while 1^1,2,3 are mostly active. 1^1,2,3 are defined in the "usual way": m\ < m\ while Irrig — mj| > 



TOg — TO2I > "^2 ~ '^1- The case > is referred to as the 'normal' mass hierarchy, while the case mg < m\ 



and — m|| > m| — m\. The case mg > m\ 

is referred to as the 'inverted' mass hierarchy. See [3, 0] for more details 
The electron neutrino (or antineutrino) survival probability is 



p — 
^ p.p. — 



E 

1=1,2,3,4 



,n / Am%L 



l-^4|C/e.|2|[/e,f Sin^ 



1<J 



(II.l) 

(II.2) 



where L is the antineutrino propagation distance (baseline), is the neutrino energy and Amf 



Current data constrain |J7e4p < few xlO ^ for lAm^J > 10 ^ eV^ 



while lAm^gl 



Ami2 ~ 7.5 X 10~^ eV^. We are interested in L < 2 km and E^ > 2 MeV so that 



(2 - 3) X 10-3 eV^ and 



< 0.1 . 



(11.3) 



Finally, |C/eiP ~ 0.7, |f/e2p ~ 0.3, and |?7e3p < 0.04. To a good approximation, for L and Ei, values of interest. 



1 - 4(1 - |f7,3p - \Up,n\Up,\' sin^ ('^^ ) - 4(1 - \Up^)\U,,\' sin^ 



iE, 



(11.4) 



where we dropped a term proportional to 



4|C/e3nt/e4p 



sm'' ^ - snr ' '-^ 



(11.5) 



which is suppressed by four powers of small mixing angles and vanishes in the limit Am^3L/4£^^ <C 1, and a term 
proportional to 



4|C/.inf/e,Psin2 ('^^ 



8 X 10 



-3 



L 



2 km 



^ '2 MeVy 



(lU,) 



This contribution is both very small for almost all baselines and energies of interest and also well-known from combined 
solar and KamLAND data. Given the intentions of this paper, its inclusion in the following discussion is of no practical 
consequence.* We have also not included matter effects, which can be safely neglected. 

In this limit, P^e is a function of two mass-squared differences {Aml^ and Am^g) and two elements of the lepton 
mixing matrix. We will parameterize |[/e3| and \Uei\ in the "standard way" (see, for example, p]): 



|f/e3|' = cos2 Ou sin' 013 



|C/e4|' =sin2 



(11.7) 



so that 



Pp. 



l-cos^ 014 sin' 2013 sin' ('^^ 

AEy 



sin' 2014 sin' f '^^ 
AE, 



(II., 



In the limit 0i4 — we recover the well-known expression for Pee at 0i3-driven reactor neutrino experiments: Pee = 
1 — sin' 2013 sin'(Ami3L/4i<^^). Note that, in spite of the fact that we have three neutrino flavors. Pee is insensitive 



To be very concrete, one can always subtract out these very small "solar" effects from the measured Pee and perform the analysis in 
the "solar-subtracted" data sample. 
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to the "sterile mass hierarchy", i.e., it cannot tell whether rrii > 1713^2.1 or 7714 < 7713.2,1. The reason is that we are not 
sensitive to the contribution Eq. pi.Sp . 

While we have information regarding A777i3, nothing is known about Arnl^. Here we concentrate on the region 
Amf^ G [I — 10] X 10~^ eV^, for different reasons. Phenomenologically, we are interested in sterile neutrino effects 
that qualitatively modify the observed oscillation pattern. A.m\^ values smaller than 10^^ eV^ start to mimick ATTijg 
effects and are not considered, while A77ii4 values much larger than 10"^ eV^ lead to averaged out effects at both the 
near and far detectors and are hence less prominent. A detailed study of the effect of LSND/Mini-BooNE- inspired 
sterile neutrinos (A77114 > 1 eV^) was recently presented in Q. Theoretically, we have in mind a simple "seesaw" 
Lagrangian for the sterile neutrinos, where mostly active neutrino masses, mostly sterile neutrino masses and active- 
sterile mixing angles are naively related via ~ C/^ (77^3/7714), where U stands for some weighted linear combination 
of the "active" mixing angles. Therefore 6*14 values around 10^^ or 10~^ are naively related to 7714 values which are 
not more than one or two orders of magnitude larger than 7713.2,1 < 0.05 eV. Experimentally, the Bugey experiment 
disfavors sin^ 26*14 values smaller than 4 x 10~^ for A77ii4 > 10^^ eV^ while cosmological considerations constrain 
sin^ 014 < 10~^ for A7nf4 > 10""'^ eV^ [1]. Both constraints are significantly alleviated for smaller values of Amf^. 
For A7ni4 — 10~^ eV^, sin^ 2^14 values as large as 10~^ are allowed by all neutrino and cosmological/astrophysical 
data [1]. 

To qualitatively understand the effect of sterile neutrinos in reactor neutrino setups we look at 



4:E^ ~ I 10-2 eV2 ) I 400 mi \ ) 



" - 1-267 ttE^ ) ( 77^ ) ( ^ ) ■ (11.9) 



For reactor antincutrinos, the oscillation length associated to Aml^ is of order the near detector distance for Aml^ 
values in the range of interest. This means that sterile neutrinos can affect Pee in the near and far detectors in distinct 
ways. Given that the sensitivity to very small values of 6*13 relies on a "near versus far" comparison, the presence of 
such sterile neutrinos can impact the reach of these experimental setups. Before proceeding with more quantitative 
results, we present a concrete qualitative example of what we mean. 

Imagine that in order to rule out a particular ^13 value we relied solely on whether the near/far ratio deviated from 
expectations. It is useful to define the observable near/far = (Afncai7A'^^oar)/(^far/A^tar): where N is the observed 
number of events and N'^ the expected number of events (in the absence of oscillations) in the near or far detectors. 

near ^ Pee (near) 
far Pee (far) 

where Pee indicates the average electron antincutrino survival probability. In the absence of oscillations near/far = 
1. Ignoring A7ni3 effects in the near detector and assuming that in the far detector the average value of 
sin^ Amj^L/'lE^ = 1/2, 

near 1 — sin^ 20i4ai4(near) 



far 1 - cos'* 6*14 sin 26'i3ai3(far) - 0.5 sin 26*14 

where is the average values of siii^ {ArnfjL/AE^) in the near or the far detector. A measurement of near/far 
consistent with one would be consistent with nonzero ^13 values satisfying 



sin^ 26*14014 (near) — cos"* 6*14 sin^ 26*13013 (far) + 0.5 sin^ 2^14 , (11-12) 
sin^ 26*14 / ai4(near) — " ^ ^ 
cos'^Ou \ oi3(far) 



• 2o/) sin^ 26*14 /ai4(near) - 0.5 \ m1Q^ 

sm 26*13 = 47— ■ (11-13 

cos* 6*14 V oi3(far) / 

6*13 effects in the far detector can be "compensated" by 6*14 effects in the near detector as long as oi4(near) is larger 
than one half and 6*14 is large enough. If ai3(far)~ ai4(near)~ 1 such an effect occurs if sin^ 2^13 ~ sin^ 2^?i4/2. In 
summary, the absence of a discrepant near/far ratio can be interpreted in one of two ways: either there are only three 
neutrinos and ^13 is very small, or there are four neutrinos and Arnl^ and ^14 are such that ^13 effects in the near 
detector and ^14 effects in the far detector "cancel!" 

The results of the above simplified analysis need to be qualified. In the case of "large" A777 i4, oi4(near) ^1/2, and 
the ambiguity is erased (this is the scenario discussed in Q). More importantly, in the case of large enough ^43 and 
9i4 (or, equivalently, once a large enough data sample is available), one expects to extract most of the information 
from distortions (or lack thereof) in the electron antineutrino energy spectrum. Ultimately, we expect that ^14 effects 
will lead to a slight loss of sensitivity to the smallest values of ^13, i.e., those close to the ^13 sensitivity boundaries 
computed in [ll,|2|. 
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Nonzero values of 614, on the other hand, may lead to potentially very unexpected results. For example, it is easy 
to see that near/far ratio may exceed one (not possible when ^14 = 0). In the case of ^13 = 0, near/far is proportional 
to 

near 1 — sin^ 2^14014 (near) 
la^ 1 -O.Ssin^ 26ii4 ' 

In this case, a near/far result different from one leads to a measurement of both sin^ 2^14 and Aml^. 



III. NUMERICAL RESULTS 



In order to study the impact of very light sterile neutrinos on 0i3-driven reactor antineutrino experiments, we 
simulate data for different values of the mixing parameters Amf3 and ^13, 14 and proceed to analyze these data 
under distinct hypotheses. The expected number of events Ni at a detector located a distance L from the source 
(assumed to be point-like) which are associated to incoming neutrino energies between Ei and Ei + AE is 



iV,(L) - NoiL) 



Ei+AE 



(j){E,)a{E,)P,e{E,)dE,, 



{111.1) 



Ei 



where 4>{E^) is the energy dependent antineutrino flux while <t{E^) is the total cross-section for Df, -\- p ^ -\- n. 
For concreteness, we use the expression for the time-averaged 4>{Ev) adopted in (see also [ll|) while we use the 
expression for a{E^) computed in [13]. Nq{L) is a normalization factor that depends on the size of the detector, the 
intensity of the source, the running time and the source-detector distance {Nq{L) oc L~^). 
We analyze our different simulated data sets by performing a simple test, where 



X^(Ami4, 6*14, Ami3, 013, a) 



^= (TV, - (1 + a)T,f 



1=1 



{my 



(iV, - (1 + a)T,f 

i=l 



far 



(III.2) 



Here Ti is the theoretically expected number of events in the i-th energy bin, given by Eq. (|nLT|l . T, depends on 

^2 /q„ jg ^YiQ error on Ni, which we assume is purely statistical and 



the oscillation parameters Am^4, ^14, Amjg, (713. 



Gaussian: 5Ni — ^/Ni. The sums are performed over all ribins energy bins in the near and far detectors, a is a nuisance 
parameter that governs how well the expected number of reactor antineutrino induced events at zero baseline can be 
predicted. 5a contains all uncertainties that are common to the expected number of events at both the near and far 
detectors, including uncertainties on the overall reactor antineutrino flux, uncertainties on the energy dependence of 
the flux, uncertainties on the antineutrino-target cross-sections, etc. No other systematic effects are considered, and, 
for the purposes of this analysis, we assume that a does not depend on the energy bin. The role of the a parameter 
is simple. Roughly speaking, when Ti/Ni — 1 is smaller than Sa, most of the statistical power in the analysis comes 
from a comparison of near detector versus far detector expectations since in the limit where either detector is "turned 
off" one can choose a such that {Ni — (1 + a)Ti)'^ is small while a'^/{Sa)'^ is order one. 

We consider two different simulated setups. In the Double-CHOOZ-like setup we choose Lnoar = 400 m, 
LfnY = 1050 m and No{L) such that one would accumulate 40,000 events in the far detector in the absence of 
oscillations. For the near detector, we set the number of unoscillated events equal to that in the far detector times 
-^farZ-^ncar-^ In the Daya-Bay-likc setup 0, we choose Lnoar = 400 m, Lfar = 1800 m and No{L) such that one would 
accumulate 75,000 events in the far detector in the absence of oscillations. As in the Double-CHOOZ-like setup, we 
set the unoscillated number of events in the near detector equal to that in the far detector times ifar/^ncar- We fix 
Sa = 3%, which is of order the overall systematic uncertainty estimated in the CHOOZ experiment il3|. Both the 
Double-CHOOZ and Daya Bay collaborations are aiming at understanding their experimental setups at a level that 
translates into, roughly, 6a ~ 1%. In what follows, we briefly comment on the impact of allowing for smaller Sa 
values. 

Two sample "data sets" for the Daya-Bay-like setup can be found in Fig. [1] It depicts the number of events 
normalized to the expected number of events in the absence of oscillations in twenty equal- width recoil positron kinetic 



t In the Double CHOOZ experiment, the near detector is expected to be smaller than the far detector. However, since most of the results 
presented here are dominated by the statistics in the far detector, we find that this discrepancy does not lead to any qualitatively distinct 
results. The same comment applies to the Daya-Bay-like setup. 



- 1.293 MeV) between 1 and 8 MeV. In both panels, sin'' 26*13 = 0.042, 
eV^. In the left-hand side, 614 = so that, as far as this observable (P, 



2.2 X 10-3 eV^, 
at the baselines and 



energy bins {Ee = 
Amj^ = 1.0 X 10 

energies of interest here) is concerned, there are no sterile neutrinos. In the right-hand side, sin^ 20i4 = 0.069. In 
the case sin^ 2^14 = 0.069 two features are noteworthy. One is that the Aml^ effects lead to visible distortions in 
the recoil electron energy spectrum both in the neat and far detectors, assuming the energy resolution is such that 
one can "see" the binning depicted in the figure. The other is that, on average, the near and far detectors point to a 
similar suppression of the expected number of events, as discussed in Sec. HH 



Daya Bay 



, ' ' I 
sin^2e^4=0.069 

^ ^1.02^ 




EglMeV) 
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FIG. 1: Number of events per energy bin in the Daya-Bay-like setup, normalized to the expected number of events in the 
absence of oscillations. Error bars are statistical only. The "data" correspond to sin^ 2613 = 0.042, Amfs = 2.2 x 10"^ eV^, 
Am^4 — 1.0 X 10"^ eV^ and sin^ 20i4 = (left-hand side) or sin^ 2Si4 = 0.069 (right-hand side). The grey [red] open circles 
(black closed circles) with smaller (larger) error bars correspond to "data" in the near (far) detector. The dotted [blue] line 
indicates the no-oscillation case. 



A. No Evidence for Oscillations 



In the absence of an oscillation signal, next-generation 0i3-driven experiments rule out regions of the sin^ 2^13 x Am^g 
plane that are currently allowed by all neutrino data. Such a result would severely impact planning for next and 
next-next generation neutrino experiments. It would, for example, reveal that the NOt^A experiment (l^ ] cannot 
determine the neutrino mass hierarchy and strengthen the case for building a muon storage ring (neutrino factory) 

In order to study the impact of sterile neutrinos, we simulate data (as described above) consistent with no oscillations 
(013 = 9i4 = 0) and analyze it under two distinct hypotheses: (i) there are no light sterile neutrinos (as far as the 
setups in question are concerned, this is equivalent to ^14 = 0) and (ii) there is a fourth neutrino mass state with 
Amf4 e [1,10] X 10"^ eV^ and sin^ 26*14 < 0.1. In either case, we restrict Ami3 £ [2,3] x 10"^ eV^, as dictated by 
current neutrino data.^ 

Figure [2] depicts the region of parameter space ruled out at the 2a confidence level in the Double-CHOOZ-like 
(left-hand side) and Daya-Bay-like (right-hand side) setup. The darker continuous boundaries are obtained under the 
hypothesis that ^14 = 0. Note that, in spite of the simplified nature of our analyses, our results agree qualitatively 
with those in [U, [I]- The lighter [red] dashed boundaries are obtained once is marginalized over the "allowed" 
sin^ 26*14 X Ami4 parameter space. 

Throughout, we perform two different "types" of data analysis. In one case we consider a simple counting experiment 
(nhins = 1), *-e-, one counts how many electron antineutrino candidate events appear in the near and far detectors 
and compares these numbers against expectations. In this case, the ability of sterile neutrinos to "mask" ^13 effects is 



t Our Am^g window agrees with the 90% confidence level allowed range quoted in the particle data book [lf| and is much wider than the 
most recent MINOS result, Am^^ e [2.17,2.69] X 10"^ eV^ at the 2cr level 
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sin^2e 



FIG. 2: Region of the sin^ 2^13 x Am^s parameter ruled out, at the 2a confidence level in the Double-CHOOZ-like setup (left- 
hand side) and the Daya-Bay like setup (right-hand side). In the top panels we depict the result of a "total rate" experiment 
('T'bins = 1) while in the bottom the "data" is subdivided into 20 recoil positron energy bins (ribins = 20). The region to the 
right of the black, continuous line is ruled out under the hypothesis that there are no sterile neutrinos (^14 = 0). The region to 
the right of the gray (red) dashed curve is ruled out under the hypothesis that Aml^ € [1, 10] x 10~^ eV^ and sin^ 26'i4 < 0.1. 



optimal. The reason is simple. For a given value of sin^ 2^13, Aml^ one can "always" find a value of sin^ 26*14, Amf^ in 
the allowed region such that the number of events in the near detector differs from expectations as much as the number 
of events in the far detector. In this case, one cannot rely on the near/far comparison to "measure" the expected 
number of neutrino events and the sensitivity is governed by external uncertainties {Sa parameter in Eq. (jIII.2p ). 

The other case under consideration is nbins — 20, i.e., we take into account not only the overall suppression of 
the electron antineutrino flux but also potential distortions of the recoil positron energy spectrum. 20 recoil positron 
energy bins between 1 and 8 MeV corresponds to AE — 350 keV. This is slightly wider than the expected energy 
resolution at both Double-CHOOZ [H and Daya Bay In this case. Fig. [5] reveals that the ability to exclude 
sin^ 2^13 values is not severely compromised by the light sterile neutrino hypothesis. The ribins = 20 case is, perhaps, 
a more faithful estimate of the results one would obtain with a realistic detector simulation. Detector and energy 
dependent systematic effects (not included in our analyses), however, tend to reduce the power of the binned analysis 
compared to overall flux one. 

Figure [3] depicts the allowed region of parameter space in the sin^ 2(?i3 x sin^ 2614 plane assuming the data in the 
Double-CHOOZ-like setup (left-hand side) and in the Daya-Bay-like setup (right-hand side) are consistent with no 
oscillations and after is marginalized over the allowed values of Am^g and Aml^ . In the counting experiment case 
(horizontal-vertical grey [red] hatching) one sees that when larger sin^ 2^14 values are considered, larger 613 values 
are consistent with no oscillations. This is in agreement with the estimate made in Eq. I|II.13|) . which can be roughly 
translated as follows. By allowing different values of Aml^, the near/far ratio, on average, can be set to one for all 
6*14 values satisfying sin^ 2^14 < 2sin^ 26*13, assuming ai3(far) ^ 1. Note that if 6*13 were set to zero, the analysis of 
a counting experiment consistent with no oscillations cannot rule out sin^ 2^14 = 0.1. This is due to the fact that 
for nbins ~ 1 and large enough Aml^ 0i4-driven oscillations average out at both the near and the far detectors. In 
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this case the sensitivity to sin^ 26*14 is dominated by Sa. Since that was set to 3%, in quaHtative agreement with 
the original CHOOZ experiment, the sensitivity to sin^ 29 is similar to the sensitivity of the CHOOZ experiment to 
averaged out oscillations, sin^ 29 ^ 0.1. For smaller values of Sa, values of sin^ 29u larger than several percent are 
ruled out even in the ribins = 1 case. 




-1 



sin^20^3 sin^20^3 



FIG. 3: Allowed region of parameter space (2a confidence level) in the sin^ 29is x sin'^ 2di4 plane assuming the data in the 
Double-CHOOZ-like setup (left-hand side) and in the Daya-Bay-like setup (right-hand side) are consistent with no oscillations. 
The (red) horizontal-vertical hatching indicates the result of a "total rate" experiment (ubins = 1) while the black criss-crossed 
hatching corresponds to the "data" subdivided into 20 recoil positron energy bins (ribins = 20). 

In the case of the more finely binned data (20 bins, criss-crossed black hatching) the absence of distortions in the 
energy spectrum at both detectors prevent large 9i3 or 9i4 values. Since 0i4-driven oscillations lead, in the near 
detector, to potentially visible distortions of the antincutrino energy spectrum even for the largest considered value of 
Am^4, the upper bound on sin^ 26'i4 is stronger than that on sin^ 2^13. For smaller values of Sa the results associated 
to nbins = 20 do not change qualitatively, while those associated to ribins ~ 1 start to approach the ribins — 20 case 
when da ~ 1%. 



B. Evidence for Oscillations ( Sia-Driven) 

If 9i3 is large (sin^ 26*43 — fewx 10^^), one expects a statistically significant disappearance of electron antineutrinos 
in the far detector of 0i3-driven reactor antineutrino experiments. In this case, one expects to not only reject the 
^13 = hypothesis but also measure ^13 (and, to a much lesser extend, Anii^). The result of such a measurement is 
also affected by whether one hypothesizes the presence of light sterile neutrinos. 

As in the previous subsection, we consider both the ribins = 1 and ribins = 20 case in order to highlight the effect of 
the light sterile neutrino hypothesis. Fig. 2] depicts the allowed region of parameter space extracted in the Daya-Bay- 
like setup assuming Amj^ = 2.2 x lO'^ eV^ sin^ 26*13 = 0.042 and 6*14 = (the "data" depicted in Fig. [H left)). The 
case ribins — 1 ("-bins = 20) is depicted in the left (right) panel. Very similar results apply for the Double-CHOOZ-like 
setup. § 

Figure m reveals that, in the case of a simple counting experiment, whether or not one allows for a light sterile state 
significantly affects the precision with which ^13 can be measured. In the ribins = 20 case, on the other hand, allowing 
for the existence of a light sterile neutrino does not significantly impact the precision with which ^13 (and Aml^) 
is measured. The reason for this is simple. In the case of 1 bin, a larger or smaller value of sin^ 2^13 can be made 
consistent with the data if it is accompanied by a large enough hypothetical sin^ 29i4 value. In the case ribins = 20, 
large values of sin^ 2^14 are ruled out by the lack of distortion in the near (and far) detector recoil positron energy 
spectrum. This phenomenon is clearly seen in Fig.[S]Jright), which depicts the allowed region of the sin^ 2^14 x sin^ 2^43 



§ Henceforth, we restrict our discussions to the Daya-Bay-like setup, keeping in mind that parallel results for the Double-CHOOZ-like 
setup are qualitatively very similar. 
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FIG. 4: Allowed region of parameter space {2a confidence level) in the Arriyj, x sin^ 2Si3 plane assuming the data in the 
Daya-Bay-like setup are consistent with Ami3 = 2.2 x fO~^ eV^, sin^ 29i3 = 0.042 and 614 = 0. The light [green] solid densely 
shaded region corresponds to analyzing the data assuming 6114 = 0, while the dark hatched region is obtaned if one allows for 
a fourth light neutrino. On the left panel a "total rate" experiment is performed (wbins = 1) while in the right panel the data 
was subdivided into 20 recoil positron energy bins (ribins = 20). 



parameter space (once one marginalizes over the allowed Am^g and Ato^4 values) . It is interesting to note that in the 
f^bins = 1 case sin^ 2^14 values as large as 0.1 are allowed at the two-sigma confidence level. The reason for this is that, 
for large enough Aml^ and ribins = 1, the sensitivity to sin'^ 2^?i4 is dominated by Sa, as discussed in the previous 
subsection. For smaller values of Sa we find that "large" values of sin^ 26*14 are ruled out even in the Tibins = 1 case. 
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FIG. 5: Allowed region of parameter space in the Ami4 x sin^ 26*14 plane (left) and the sin^ 2^14 x sin^ 2^13 plane (right) 
assuming the data in the Daya-Bay-like setup are consistent with Amig = 2.2 x 10~^ eV^, sin'^ 2^13 = 0.042 and 6114 = at the 
2a confidence level. On the right panel the (red) horizontal-vertical hatching indicates the result of a "total rate" experiment 
('T'bins = 1) while the black criss-crossed hatching corresponds to the "data" subdivided into 20 recoil positron energy bins 
(^^bms = 20). On the left panel nbins = 20. 



In the case nbins = 20, the absence of a distorted positron energy spectrum in the near detector rules out sin^ 2^14 
values larger than 10~^ or so. Fig. [Sl^left) depicts the region of the Am,^^ x sin^ 2^14 plane ruled out at the two sigma 
level assuming the data in the Daya-Bay-like setup are consistent with Am^g = 2.2 x 10"'^ eV^, sin^ 26*13 — 0.042 and 
^14 = 0. Such a result would improve the current constraints on ^14 by about an order of magnitude for the values of 
Ami4 highlighted here. 
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C. Evidence for Oscillations (S13 and Si4-Driven) 



If 9i4 is nonzero (and large enough) and /^m\^ ^ few x 10^^ cV^, 6'14-efFects will produce non-trivial, distinct 
effects in the near and far detectors, as depicted in Fig. fright). In this case, allowing for the presence of a light 
sterile neutrino during the data analysis would prove to be more than a choice — it would be necessary in order to 
obtain a proper fit to the data of the reactor neutrino experiments under investigation here. 

Figure [5] depicts two two-parameter 2a allowed regions of parameter space if the data in the Daya-Bay-like setup 
were consistent with Am^g = 2.2 x lO'^ eV^ sin^ 26*13 = 0.042, Am?4 = 0.01 eV^ and sin^ 26*14 = 0.03. Figure E^left) 
depicts how well sin^ 2^13 (and Am^3) can be measured in the case nbins = 1 and ribins = 20 after one marginalizes 
over sin^ 26*14 and /S.m\^. These should be compared to Fig. [4l It is clear that in the ribins = 1 case the fact that 
sin^ 26^14 is nonzero renders a "rates-only" measurement of sin^ 2^13 harder. In the nbms = 20 case, on the other hand, 
the impact of a non-zero sin^ 2^14 when it comes to measuring sin^ 2^13 is small. The reason for this is that in the 
'T-bins = 20 case distortions in the near (and, to a lesser extent, in the far) detector determine Arnf^ and sin^20i4, 
allowing the near/far comparison to "fully contribute" to the measurement of ^13. 
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FIG. 6: Allowed region of parameter space in the Aml^ x sin'^ 2^13 plane (left) and the sin^ 2614 x sin'^ 29is plane (right) 
assuming the data in the Daya-Bay-like setup are consistent with Amf3 = 2.2 x 10"'^ eV'^, sin^ 2613 = 0.042, Ami4 = 0.01 eV^, 
and sin^ 2^14 = 0.03 at the 2a confidence level. On the left panel, the light [green] solid densely shaded region corresponds to 
fibins = 20, while the dark hatched region corresponds to ribins = 1. On the right panel, the (red) horizontal-vertical hatching 
indicates the result of a "total rate" experiment (ribins = 1) while the black criss-crossed hatching corresponds to the "data" 
subdivided into 20 recoil positron energy bins (ribins ~ 20). 



Figure [BKright) depicts how well sin^ 2^13 and sin^ 2^14 can be measured in the case ribins — 1 and ribins = 20 
after one marginalizes over both mass-squared differences. In the case ribins ~ 1, virtually no constraint can be set 
on sin'^ 26*13. Curiously enough, if sterile effects were not included (6*14 = 0), one would be able to establish that 
sin^ 2^13 7^ 0. In this case, however, the wrong hypothesis in the data analysis would point to an allowed range for 
sin^ 2^13 that is slightly less than its real value (this is true at around the la level). For larger "true" values of ^14 
this effect is more pronounced. As discussed earlier, in the case ribins — 20 one is able to obtain a precise measurement 
of both mixing angles. 

For smaller values of true sin^ 29i4, the impact of sterile neutrinos is, of course, less pronounced. Figure [7] depicts 
two two-parameter 2a allowed regions of parameter space if the data in the Daya-Bay-like setup were consistent with 
Am^3 = 2.2 X 10-3 eV^ sin^ 26*13 = 0.042, Am^4 = 0.1 eV^, and sin^ 26*14 = 0.0091. In this case, even in the ribins = 1 
case one can establish that sin^ 26*13 7^ (at the 2a level). In the ribins = 20 case, a reasonable measurement of both 
sin^ 26*13 and sin^ 26*14 can be performed — at the 2a level, sin^ 26*14 e [0.002, 0.015] and sin^ 26*13 S [0.03, 0.055], and 
correlations are small. 

We conclude this section with a short comment on the case 6*13 = and 614 7^ large. As far as measuring 
sin^ 2^14 is concerned, the situation here is qualitatively similar to the previous two cases displayed above. Due to 
the cancellation effects discussed earlier, an ribins — 1 analysis fails to severely constrain either ^?i3 (similar to the 
situation depicted in Fig.l^^top)) and ^14. An ribins — 20 analysis, on the other hand, allows one to not only measure 
sin^ 2^14 and Arnf^ with good precision but also constrain sin'^ 26*13 significantly, as in Fig. [^Kbottom). 
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FIG. 7: Allowed region of parameter space in the Ami3 x sin^ 26^13 plane (left) and the sin^ 2di4 x sin^ 26i3 plane (right) 
assuming the data in the Daya-Bay-like setup are consistent with Amig — 2.2 x 10"'^ eV^, sin^ 26i3 — 0.042, Anii^ — 0.1 eV^, 
and sin'^ 29i4 — 0.0091 at the 2a confidence level. On the left panel, the light [green] solid densely shaded region corresponds to 
'T-bins = 20, while the dark hatched region corresponds to ribins = 1. On the right panel, the (red) horizontal-vertical hatching 
indicates the result of a "total rate" experiment (nbins = 1) while the black criss-crossed hatching corresponds to the "data" 
subdivided into 20 recoil positron energy bins (ribins ~ 20). 

IV. DISCUSSION AND CONCLUSIONS 

We have discussed the impact of very light sterile neutrinos on the interpretation of next-generation 0i3-driven 
reactor antineutrino experiments. For Aml^ values between 1 and 10 xlO~^ eV^, 0i4-driven oscillations affect data 
in the far and near detectors differently, leading to several potentially interesting effects. In the absence of a positive 
oscillation signal, the hypothesis 614 7^ negatively impacts one's ability to rule out very small 613 values as would- 
be 014-driven effects in the near and far detectors can mask 0i3-driven effects in the far detector. Similarly, if a 
positive 6'13-driven signal is observed, would-be 0i4-driven effects negatively impact ones ability to measure the value 
of sin^20i3. The above results are more or less pronounced depending on one's abihty to "see" distortions of the 
recoil positron energy spectrum. A simple counting experiment is very susceptible to hidden 0i4-driven effects while 
a binned data analysis (at both the near and far detectors) seems able to disentangle ^13 from 0i4-driven effects. 

If there is indeed a light sterile neutrino associated to Aml^ E [1, 10] x 10^^ eV^, 0i3-driven reactor antineutrino 
experiments should be able to see large spectral distortions in the near detector, in which case one can determine 
not only 614 but also Am^4. It is curious to note that, for Aml^ values close to 0.1 eV^, the roles of the near and 
far detectors as far as studying 0i4-driven effects are reversed compared to those associated to studying 0i3-driven 
effects. For these mass-squared differences, oscillation effects average out in the far detector (even in a finely-binned 
analysis). This allows one to determine the "averaged-out" neutrino flux and hence helps extract the value of ^14 
from the depths of the oscillation minima. In the absence of 0i4-driven effects, we estimate that values of sin 2^r4 as 
small as 1 percent can be ruled out with a binned analysis (Fig. [5l[left)). 

Sterile neutrino effects at 0i3-driven reactor experiments have been considered in the past, but in a different mass- 
squared difference regime. The authors of foj recently discussed the effect of sterile neutrinos related to a potential 
3 + 2 solution to the LSND anomaly.^ There, the new mixing angles and mass-squared differences lead to averaged 
out effects at both the near and far detectors. Here, on the other hand, we concentrate on a lower mass-squared 
difference regime where this is not the case. 

If there are indeed such light sterile neutrinos, it is likely that the most sensitive terrestrial probes of their existence 
are the setups discussed here. This is not a coincidence: we are concentrating on values of Am\4 where oscillation 
effects in the near detector are optimal. Other near-future probes include long-baseline accelerator-based neutrino 
oscillation experiments. However, Amj4-driven effects at next-generation 0i3-driven appearance long-baseline exper- 



^ After MiniBooNE data became available, it was pointed out that a 3+2 solution to all short baseline appearance data exists. This 
solution, however, is disfavored by disappearance searches sensitive to the same oscillation frequencies |18| . 
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iments are proportional to the product of the squares of two potentially small mixing angles (P^e c>c 91^024^, using 
the notation of Furthermore, next-generation studies of muon neutrino disappearance at similar setups are also 
unlikely to achieve sensitivity competitive with the one discussed here. Note that Amf4 is large enough that it leads 
to averaged-out effects at the far detector but too small to mediate observable effects in the near detectors of next- 
generation long-baseline experiments.** Detailed recent analyses of next-generation oscillation probes of light sterile 
neutrinos, concentrating on Aml^ values above 0.1 eV^, can be found in |19l. [20l. [2ll. 12^ . 

It is important to appreciate that this range of sterile neutrino parameters is not motivated by any existing data, nor 
does it help address any existing outstanding issue in fundamental physics. On the other hand, as already emphasized, 
little is known about sterile neutrinos. Evidence for sterile neutrinos at any mass range would qualitatively change 
our understanding of particle physics (and probably shed light on the mechanism behind light neutrino masses). 

Finally, light sterile neutrinos qualify as an example of non-standard physics that may appear in next-generation 
reactor neutrino experiments. They can not only lead to nontrivial oscillation patterns in the data, but also modify the 
interpretation of reactor antineutrino data when it comes to measuring or constraining 9-\^. Our results also highlight 
the fact that new and interesting results may come out of the data in the near detector [23| , which may provide more 
information than the measurement of the "L = 0" reactor antineutrino flux. We conclude by re-emphasizing that our 
simulations and analyses are not aimed at realistically describing experimental setups or quantitatively gauging their 
reach. Qualitatively, however, our results capture the nontrivial impact of light sterile neutrinos at 0i3-driven reactor 
antineutrino experiments. We hope that our findings will prompt the collaborations to pursue quantitative estimates 
of the impact of the very light sterile neutrinos introduced here. 
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